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ABSTRACT 

Here we discuss the mechanical feedback that massive stellar clusters provide 
to the interstellar medium of their host galaxy. We apply an analytic theory 
developed in a previous study for M82-A1 to a sample of 10 clusters located in 
the central zone of the starburst galaxy M82, all surrounded by compact and 
dense HII regions. We claim that the only way that such HII regions can survive 
around the selected clusters, is if they are embedded into a high pressure ISM and 
if the majority of their mechanical energy is lost within the star cluster volume via 
strong radiative cooling. The latter implies that these clusters have a low heating 
efficiency, 1], and evolve in the bimodal hydrodynamic regime. In this regime 
the shock-heated plasma in the central zones of a cluster becomes thermally 
unstable, loses its pressure and is accumulated there, whereas the matter injected 
by supernovae and stellar winds outside of this volume forms a high velocity 
outflow - the star cluster wind. We calculated the heating efficiency for each of 
the selected clusters and found that in all cases it does not exceed 10% . Such low 
heating efficiency values imply a low mechanical energy output and the impact 
that the selected clusters provide to the ISM of M82 is thus much smaller than 
what one would expect using stellar cluster synthetic models. 

Subject headings: galaxies: individual (M82) — galaxies: star clusters — HII 
regions — ISM: bubbles — kinematics and dynamics 

1. Introduction 

It is a common believe that massive star clusters return a significant fraction of their 
stellar mass to the interstellar medium (ISM). This is thought to be done in a violent manner 
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that deeply affects tfie structure of the interstellar gas and in the case of starburst galaxies, 
this may even result in the channeling of the processed material into the intergalactic space 
(see, for instance, Tenorio-Tagle et al. 2003, Cooper et al. 2008, and references therein). 
The general concensus is that within the volume occupied by superstar clusters (SSCs), 
the kinetic energy supplied by massive stars in the form of stellar winds and supcrnovac 
explosions is there in situ thermalized. This results into a high temperature (T ~ lO^K) 
plasma, with a large thermal pressure that highly exceeds that in the ambient ISM, and this 
provokes the exit of the thermalized ejecta out of the cluster as a supersonic star cluster 
wind (Chevalier & Clegg, 1985). The cluster winds shape the ISM by generating large-scale 
superbubbles. These shock and displace the surrounding ISM while locking it into large 
expanding shells able to cool down by radiation in a short characteristic time scale, while 
the much lower density shock-heated wind gas, which fills the supcrbubblc interior, remains 
hot for a considerably longer time (Weaver et al. 1977; Mac Low & McCray, 1988; Tenorio- 
Tagle et al. 2006) and promotes the growth of the superbubble. This shocked wind plasma 
has been detected around OB-associations and stellar clusters as a soft X-ray emitter (sec, 
for example, Chu et al. 1995; Stevens & Hatwell 2003, Sihch et al. 2005 and references 
therein), whereas the outer shells have been traced in 21 cm (Puche et al. 1992, Ehlerova et 
al. 2004) or as photoionized engulfing filaments, if in presence of a strong Lyman continuum 
radiation (Meaburn, 1980; Lozinskaya, 1992 and references therein). The size and interior 
pressure of superbubbles in the case of an homogeneous interstellar gas distribution (see Mac 
Low & McCray, 1988; Bisnovatyi-Kogan & Sihch 1995 and references therein) are: 
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where 7 = 5/3 is the ratio of specific heats, Rgb is the outer shell radius, Psb is the pressure 
in the shock-heated wind region. Lout is the star cluster mechanical energy output, pisM is 
the interstellar gas density, and t the evolutionary time. The superbubbles are supposed to 
expand until they acquire pressure equilibrium with the surrounding medium {Pgh = Pism), 
when: 
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where rijsM is the interstellar gas number density, L^s is the mechanical energy output in 
units of 10^^ erg s~^ and Oio is the sound speed in the interstellar medium in units of 10 km 



However, as noticed by Sihch et al. (2007), this cannot be the whole story. The observed 
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properties of the HII region associated to the massive (1.3x10^ M©), young (age ~ 6 Myr) 
and thus powerful (with a mechanical luminosity, Lmech ~ 2.5 x 10^° erg s~^) super star 
cluster M82-A1 (Smith et al. 2006) are not consistent with the interstellar bubble model 
(equations [T] - [3]). The associated minute, low mass {Mhii ~ 5 x 10^ Mq), although dense 
{hhii ~ 1800 cm~^), HII region presents a radius {Rhii ~ 4.5 pc) much smaller than that 
predicted by equations ([Hand [3]). It seems surprising that M82-A1 and other young and 
massive clusters in M82 (Melo et al. 2005), NGC 3351 (Hagele et al. 2007) and in other 
galaxies are surrounded by compact, low mass (see Table 2 below) HII regions despite the 
powerful mechanical energy output predicted for the clusters by stellar evolution synthesis 
models. 

Silich et al. (2007) suggested that in this case clearly only a fraction of the star cluster 
mechanical luminosity is converted into the energy of the outflowing plasma whereas the 
rest ought to be lost due to strong radiative cooling. They also developed an analytic and 
semi-analytic model, which led to obtain the value of the heating efficiency rj, the parameter 
which links the star cluster mechanical luminosity with the actual thermal energy that is 
deposited into the star cluster volume. The models reveal the value of the heating efficiency 
by fitting the ionized gas number density and radius of the compact HII region detected 
around a massive star cluster. The results led to a low heating efficiency [rj < 10%) in the 
case of M82-A1. 

Here we extend the analysis of Silich et al. (2007) to a sample of 10 SSCs selected from 
the list of Melo et al. (2005), in order to reveal their heating efficiency, and thus the energy 
that these clusters return to the ISM of M82. 

Section 2 and 3 discuss the hydrodynamics of SSCs, the model assumptions and present 
the equations used in order to obtain the heating efficiency. The sample of the selected 
clusters is presented in section 4. We apply our model to each of the selected clusters and 
discuss our results in section 5. 



2. The heating efficiency in SSCs 

The hydrodynamics of the matter returned by stellar winds and supernova explosions 
within the cluster volume has been approximated assuming first that the sources are equally 
spaced within a spherical volume of radius Rsc- In the pioneer adiabatic approach of Cheva- 
lier & Clegg (1985), the kinetic energy supplied by the evolving massive stars, Lmech, has 
been assumed to be completely converted into thermal energy of the hot plasma. The strong 
pressure gradient generated by the deposited matter, forces then the gas velocity to increase 
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almost linearly from km at the star cluster center to its sound speed at the cluster edge. 
Once the gas streams out of the cluster, it rapidly acquires its terminal speed {vaoo ~ '^csc), 
while its density and temperature drop as r~^, and r~^^^, respectively. 

More recently, Sihch et al. (2004), Tenorio-Tagle et al. (2005, 2007) and Wiinsch et 
al. (2008) recognized that the adiabatic assumption is not vahd in the case of massive and 
compact star clusters and developed a radiative star cluster wind model, which takes into 
consideration the energy losses that occur in the hot thermalized plasma. They found a 
threshold line, Lcrit{Rsc) in the L^ng^h — Rsc parameter space. The radiative solution is 
in excellent agreement with Chevalier & Clegg's results in the case of low mass clusters, 
when Lmech Lcrit- However, strong radiative cooling modifies essentially the temperature 
distribution outside of the cluster when the star cluster mechanical luminosity approaches the 
threshold value, L^ech < -^crit- When the mechanical luminosity of the considered clusters 
exceeds the threshold value, Lcrit, catastrophic cooling sets in within the central zones of 
the cluster, what results into a bimodal flow regime. In this case the stagnation radius, Rgt, 
moves out of the cluster center and splits the cluster volume into two distinct zones. In the 
inner zone, r < Rgt, strong radiative cooling promotes frequent thermal instabilities in the 
injected gas, reducing significantly the pressure gradient and thus the outward acceleration. 
Strong radiative cooling thus leads to the accumulation of the matter injected within the 
volume defined by the stagnation surface. In the outer zone, Rgt < r < Rsc, despite radiative 
cooling, the pressure gradient remains sufficient to drive the injected matter away from the 
cluster, as a strongly radiative stationary wind. 

The detail physics during the thermalization process are however not well understood. 
In the original paper of Chevalier & Clegg (1985) it was assumed that the amount of the 
deposited thermal energy per unit volume, Qg, is identical to the rate of mechanical energy 
released by massive stars: Qe — Qmech- However Stevens & Hartwell (2003) found that this 
assumption is not in good agreement with the spectra of the diffuse X-ray emission detected 
in a number of nearby massive clusters. Bradamante et al. (1998), and Recchi et al. (2001), 
who studied the chemical and dynamical evolution of blue compact galaxies, also claimed 
that only a few per cent of the energy deposited by supernovae type II provides the energetics 
of the host galaxy ISM, while the rest is radiated away. It is therefore highly desirable to 
fink the value of the heating efficiency with stellar clusters observable quantities. 

A firm evidence for an incomplete transformation of the star cluster mechanical lumi- 
nosity into the energy of the star cluster wind was obtained by Smith et al. (2006), who 
provided detailed photometric and spectral analysis of the massive, young SSC M82-A1 and 
its associated HII region. This led them, as well as to Sihch et al. (2007), to claim that 
the energy qe represents only a small fraction of the mechanical energy provided by massive 
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stars: Qe — rjqmech, with 77 << 1. The physical justification for this parameter comes from 
the fact that strong radiative coohng may take place during the process of thermalization 
either because of an enhanced gas metallicity, resultant from SN explosions, or because of 
the large densities within the shock-heated zones between neighboring massive stars, before 
the newly injected matter joins the flow (Wiinsch et al. 2007; Silich et al. 2007). In this 
case only a fraction of the mechanical energy supplied by the collection of massive stars is 
shared by the matter within the cluster volume and thus the actual thermal energy given 
to the injected gas is smaller than that provided by the collection of massive stars. This is 
particularly important in the case of massive {Msc > 10^ Mq) and compact {Rsc ~ a few 
parsecs) clusters which present a large massive star number density, N^, and a small mean 
separation between them: AR — <^ 1 pc. 

Indeed, the multiple interactions expected between supersonic stellar winds and the 
supernovae ejecta in such compact and massive clusters are similar to those occurring in 
colliding wind binaries that lead to a shock-heated plasma, that effectively radiates in the soft 
X-ray regime. Luo et al. (1990); Stevens et al. (1992) found that in the case of colliding wind 
binaries the amount of energy radiated away from the shock-heated zone, Liost, depends on 
the binary separation. It scales as Liost ~ AR~^ in the quasi-adiabatic regime and increases 
when radiative cooling in the shock-heated zone is taken into consideration. In the case of a 
dense stellar cluster the kinetic energy placed by massive stars interacts with that deposited 
by multiple nearby neighbors. This suggests that the energy, which actually drives the star 
cluster outflow, is smaller than the total provided by the massive stars within the cluster 
volume, particularly if one accounts for the large metallicities expected from supernovae. 

In the semi-analytic models all uncertainties dealing with the distribution of massive 
stars and the collisions between nearby supersonic flows and thus the sudden loss of energy, 
are accounted for by the parameter 77, known as the heating efficiency. The fraction of energy 
that a star cluster returns to the ambient interstellar gas strongly depends on this parameter. 
Thus rj defines the mechanical feedback that star clusters provide to the ISM of their host 
galaxy. 

3. The pressure confined wind model 

The pressure confined wind model (Sihch et al. 2007) suggests that the combination 
of two factors is crucial in order to produce the compact and dense HII regions able to 
survive around powerful young clusters. These are a low heating efficiency and a large 
thermal pressure, Pism, in the surrounding ISM, what leads to a pressure confined bubble 
configuration. Thus in this model the size of the standing HII region depends critically 
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on the balance between Pjsm and the wind ram pressure at the reverse shock position 
{Pram = Pism)- this casc the structure of the outflow can be derived analytically from 
a set of equations that consider: conservation of mass, photoionization balance, pressure 
equilibrium and the fast radiative cooling that occurs within the star cluster volume and on 
its wind. Our set of equations is such that if the parameters of the driving stellar cluster: 
its mass (Msc), radius (Rsc), the number of ionizing photons {N^'~"') are known, one can 
match the model predicted radius {Rrii) and gas number density [riHjj) of the associated 
HII region with the observed values. In this approach, one can find the value of the heating 
efficiency from a nonlinear algebraic equation which relates rj with the host cluster and the 
associated HII region parameters (see Silich et al. 2007): 
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where Vaqo = C^Lsc/ MscY^'^ is the adiabatic wind terminal speed, Lsc and Msc are the 
star cluster mechanical luminosity and the mass input rate, respectively, is the value of 
the cooling function at the reverse shock radius, j3 is the recombination coefficient to all but 
the ground level, fx = 0.3 is a fiducial coefficient and ft is the fraction of ionizing photons, 
which reaches the outer standing gaseous shell. Note that equation (jl]) is only valid in the 
bimodal parameter space, i.e. it can only be applied to clusters with a mechanical power 
that exceeds the threshold value Lcru, and carries a strong implicit dependence on rj via 
the threshold mechanical luminosity, Lcrit and the star cluster wind terminal speed, Voo (see 
Wiinsch et al. 2007; Sihch et al. 2007): 
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where a = 0.28 and fii = lArriH/H is the mean mass per ion. Agt and Cst are the values of 
the cooling function and the speed of sound at the stagnation point, both are functions of 
temperature at the stagnation radius, Tst, which strongly depends also on rj. We obtain the 
value of Tst from the condition that the stagnation pressure reaches the maximum possible 
value and thus dPgt/dTst = (Tenorio-Tagle et al. 2007): 
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Note that the calculated heating efficiency does not depend significantly on the parameter 
ft- Hereafter we shall assume that ft = 0.5 and Vaoo = 1000 km s~^. 
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In order to relate the star cluster mass with the star cluster mechanical luminosity we 
use a relation, which approximates the results of Starburst 99 synthesis model for coeval 
clusters with a Salpeter initial mass function with sources between 1 Mq and 100 Mq and 
ages in the range ~ 4 Myr - 12 Myr (Leitherer et al. 1999): 



Equation (jlj) presents only a weak dependence on Lsc and thus deviations of the star clus- 
ter mechanical luminosity from the assumed constant value do not affect the final results 
significantly. 

Thus, in this approach one can obtain the heating efficiency rj directly from the observed 
parameters of the stellar cluster and its associated HII region: Msc, Rsc, N^'-^ , Rhii and 
nHii by solving equation (j4]). 



In order to learn how efficient the conversion of the star cluster mechanical luminosity 
into the wind driving energy is, one needs a sample of clusters whose masses, sizes and Lyman 
continuum radiation are known together with the radius and density of their adjacent HII 
regions. Most of these parameters can be obtained from the photometric sample of Melo et 
al. (2005) who cataloged 197 young super stellar clusters in the central zone of the galaxy 
M82. The only parameter, which is required by equation (jlj) and which Melo et al. (2005) 
did not obtain, is the density of the ionized gas in the HII regions. We obtain this quantity 
from Potsdam Multi Aperture Spectrophotometer, PMAS (Roth et al. 2005) observations at 
the 3.5 m telescope in Calar Alto. PMAS is a very versatile instrument, with several working 
modes. Here, we used its lens array (LARR) which is made out of 16 x 16 square elements. 
We observed the nuclear region of M82 using two continuous fields with the spatial sampling 
of 0'.'5 X 0'.'5 and thus, covering a field of view of 8'.'0 x 8'.'0 per pointing. Two different sets of 
observations were provided. A set of data with low spectral resolution, covering the whole 
optical spectral range was obtained in service mode on the 3rd and 4th of June 2005 while 
that with a high spectral resolution on the 2nd of February 2005. Both sets of data were 
taken under non-photometric conditions. Seeing ranged typically between 1'.'3 and 1'.'6. Line 
profiles were fitted using gaussian functions. This procedure was done in an automatic way 
using the IDL based routine MPFITEXPR implemented by Markwardt0 checking each fit. 
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4. A sample of clusters in M82 
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afterwards. A single gaussian fit was enough in order to reproduce the observed profiles. For 
each set of lines, wavelenght differences between them were fixed and the same line width 
was assumed. Then the intensity maps of and 5*// lines were produced. 

Our field of view includes region A (O'Connell & Mangano, 1978) of the central zone 
as well as a highly extincted heart-shaped region towards the west. In order to localize the 
selected clusters in the PMAS map, the resolution of the HST image was degraded to 0'.'506 
/ pixel, which is almost identical to that of PMAS (O'.'S / pixel) and a new, low resolution 
HST Ha map of M82 nuclear region was generated. As we did not have absolute astrometry 
and in order to use the two observing data sets, the low-resolution HST and the PMAS Ha 
maps were compared and displaced until reaching the highest cross correlation coefficient 
(Russ, 2002). Then the electron density, uhhi was derived (see Table 2) from a map of 
the [Sll]A6717 / [Sll]A6731 line ratio using the task temden, based on the fivel program 
(Shaw & Dufour, 1995), included in the IRAF package nebular and assuming an electronic 
temperature of 10000 K (see Figure 1). 




10 20 30 

pixel 



Fig. 1. — The ionized gas density distribution derived from the PMAS observations. The 
range of densities is coded by color scale on the right of the panel. Positions of different 
points in the field of view of PMAS are given in pixels. The position of the M82-A1 cluster 
is marked with a cross symbol. 
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We use the electron density values and the photoionization balance equation in order 
to estimate the masses of the associated HII regions: 

Mhh = (9) 

where /ij is the mean mass per ion and (3 = 2.59 x 10~^^ cm~^ s~^ is the recombination 
coefficient to all but the ground level. 

84 out of the 197 clusters found by Melo et al. (2005) are located in the area observed 
with PMAS. From the large sample of young SSCs cataloged by Melo ct al. (2005) we 
have selected a subsamplc which follows the criteria that the radius of the Hll region (the 
one defined in the Ha HST images) lies clearly outside the volume ocupied by the SSCs 
themselves (radius taken from the continuum HST images) . In this way we selected a total 
of 21 objects. We then compared our list of cluster candidates with that of Mayya et al. 
(2008), who used the HST Advanced Camera for Surveys (ACS) and selected only those 
sources, which were simultaneously detected in three different (B, V and I) filters. Only 10 
counterparts for our 21 candidate clusters were found in the list of Mayya et al. (2008). 
We have selected these as genuine clusters for our further discussion. Figure 2 presents the 
location of the selected clusters within the galaxy and also outlines the area in the central 
zone of M82, which was observed with PMAS. 

Table 1 presents the identification of the selected clusters. Here the first column marks 
the clusters in our list, columns 2 and 3 provide the star cluster and the M82 zone identifi- 
cation in the sample of Melo et al. (2005) and column 4 lists the identification number in 
the sample of Mayya et al. (2008). 

Table 2 presents the star cluster masses and radii (columns 2 and 3, respectevely) and 
the number of Lyman continuum photons (column 4) taken from Melo ct al. (2005). The 
number densities, radii and masses of the associated Hll regions arc given in columns 5, 6 and 
7, respectively. Column 8 presents the calculated values of the heating efficiency and column 
9 - the output mechanical luminosity normalized to the star cluster mechanical luminosity, 
Lsc, predicted by the Starburst 99 synthetic model. 

Note that the cluster radii fall into a narrow size interval, 2 < Rsc < 6 pc whereas their 
masses vary from 2 x 10^ M0 to 8 x 10^ M©. In all selected cases the resulting masses of 
the associated HII regions do not exceed a few thousand solar masses, just as in the case of 
M82-A1 whose stellar mass is ~ 10^ M© and its associated HII region has only ~ 5000 M© 
(Smith et al. 2006). 
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Fig. 2.— The HST WFPC2 image of the central zone of M82 taken in the F656N filter. 
The selected clusters are shown as green dots. The position of the M82-A1 cluster is used 
as a reference point and is marked with a red plus. Dashed lines outline the field of view of 
PMAS. 
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Table 1: Selected clusters 
Cluster ID zone ID 

Melo et al. (2005) Mayya et al. (2008) 



(1) 


(2) 


(3) 


(4) 


1 


20 


sw 


178N 


2 


52 


SE 


47N 


3 


12 


SE 


13N 


4 


58 


SE 


34N 


5 


14 


SE 


20N 


6 


59 


SE 


45N 


7 


3 


NE 


71N 


8 


72 


SE 


200N 


9 


81 


SE 


32N 


10 


75 


SE 


28N 



Table 2: Parameters of the selected clusters and their HII regions 



Cluster 


Msc 


Rsc 


^SC 




^HII 


Mhii 


V 






(10-> Mq) 


(pc) 


(104« s-1) 


(cm^"') 


(pc) 


(Mq) 


% 


% 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1 


0.35±oi6 


4.03 




7g9±76 


5.64 


381.6 


7.8±2-i 


317±i.06 


2 


0.4o±o.i6 


4.03 


34 0±2i.o 


950^^1 


4.83 


1479.3 


7.0±i-9 


2.46±o.93 


3 




4.83 


17.7±9-^ 


706^^1 


5.64 


1036.2 


5.0±ii 


0.83±0-34 


4 




3.22 


19.0±7-9 


953±8i 


4.03 


824.1 


5.3±i-2 


104±0.44 


5 


1.30±i-0 


3.22 


15.2±8-4 




5.64 


944.8 


5.2±i-2 


0.7i±o.3i 


6 


4.00±3-7 


4.03 


24.o±i2.o 


886=^1^5 


4.83 


1119.6 


4_0±o.8 


0.23±oi2 


7 


2.19±0.47 


3.22 


58.0±22.o 




4.83 


3109.4 


4_3±o.8 


0.33±oi5 


8 


145±o.33 


2.42 


16.6=^^-^ 


1146±76 


3.22 


598.7 


43±o.9 


034±o.i8 


9 


3.60±2i 


3.22 


18.0±8-2 


850^145 


4.03 


875.3 


3_7±o.6 


0.i7±o.08 


10 


2.4o±2.2 


2.42 


24.o±i5.o 


1163±82 


4.83 


853.0 


5.6±i-7 


0.53±0-28 



• Parameters of the clusters (columns 2, 3 and 4), radii of the associated HII regions 
(column 6) and uncertainties in their determination are taken from Melo et al. (2005). 

• Ionized gas density (column 5) was derived from PMAS observations. 

• One pixel uncertainty (±0.81 pc) was adopted in the determination of all radii. 
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5. Results and discussion 

Each of the selected clusters (see Table 2) is surrounded by a compact HII region and 
has all attributes required by our model. We solve Equation (j4]) by iteration with the relative 
accuracy ^rj/rj < 10~^. Equations ([5]), ([6]) and ([7]) were used every time when the iteration 
procedure requires new values for the threshold luminosity, Lent, cooling function, A^^, and 
the star cluster wind terminal speed, Voo. Our results for each of the considered clusters are 
shown in Table 2 (column 9). We use the error propagation equation (Bevington & Robinson, 
2003) in order to calculate the errors provided by the uncertainties in the determination of 
the input parameters of the model: Msc, Rsci N^'^ , Rhii and uhh- Unfortunately, the 
uncertainties in the determination of the star cluster radii and sizes of the HII regions are 
not presented in the original paper of Melo et al. (2005). We take a 1 pixel (±0.81 pc) as a 
conservative estimate for the uncertainties in the measured radii. 

The results of the calculations are presented in Figure 3, where the heating efficiency 
is presented as a function of star cluster radii and masses. It seems that there is a trend on 
panel b for the heating efficiency to be larger for less massive clusters. However this must 
be confirmed with better sets of input data. We also suggest for a future analysis that the 
star cluster stellar density may be a better input parameter, which combines the two major 
observables, the star cluster mass and radius, into a single parameter. 

0,101 It _ 

0,09 - 

0.08 I y. I - 

0.07 '- I a 1 - 

0-06 - T - 

~ 0,05 - T n 1 - 

0,04- ' ' 3 i 1 
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O.oa - 
0.01 - 

n no F - 

I 2 ^, A f] 6 

Rsc(pc) 

Fig. 3. — The calculated heating efficiency. Panels a and b present the star clusters heating 
efficiency as a function of star cluster radius and mass, respectively. 

Figure 3 shows that the heating efficiency does not exceed 10% for all clusters in our 
sample. This implies that our massive and compact stellar clusters have a much reduced 
outflow velocity and negative feedback into the ambient ISM than what one would expect 
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using synthetic models. Indeed, the mechanical energy output rate is: 

Lout = \MoutVL (10) 

where the star cluster wind terminal speed, Voo, is defined by equations ([7]) and ([6]) and the 
mass output rate, Mout, is ( Wiinsch et al. 2007): 

Mout = Msc(^\" . (11) 

The critical luminosity, Lcrit, is defined by equation (111). The fraction of the injected me- 
chanical energy. Lout/ Lsc : that a cluster returns to the ambient ISM thus is: 

Lout _ f Lcrit\ ^ f KxD 

Lsc V Lsc J \ Vaoo 

The calculated mechanical energy output does not exceed a few per cent of the star 
cluster mechanical luminosity for all selected clusters (see Table 2). Only in this way the 
shock-heated matter driven out as a cluster wind can cool rapidly to T < 10'^ K and be 
photoionized while a high ambient pressure prevents its expansion into the surrounding 
interstellar medium. 

The implication of our results, when compared with the recently inferred (Strickland 
& Heckman, 2009) net efficiency of supernova and stellar wind feedback in the nucleous of 
M82 (> 30%), is that there is a phase, a time during which massive and compact clusters 
have a low heating efficiency and undergo a bimodal hydrodynamic solution returning to 
the ISM of their host galaxy only a small fraction of mass and mechanical energy released 
inside the star cluster volume. Here we suggest that the selected young, massive clusters pass 
through such special phase in their hydrodynamical evolution, highlighted observationally 
by the presence of a compact HII region. Indeed, the relevant cluster parameters such as: 
the energy and mass deposition rates, the mean separation between nearby energy sources 
and the chemical composition of the injected matter - all change with time. This must lead 
to important changes in rj and thus to large displacements of the threshold luminosity and 
noticeable changes in the rates of mass, Mout, and energy. Lout, which a star cluster returns 
to the ISM. The time evolution of t] will be the subject of a forthcoming communication. 
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